Conductance quantization in graphene nanoconstrictions with mesoscopically smooth 

but atomically stepped boundaries 
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We present the results of million atom electronic quantum transport calculations for graphene 
nanoconstrictions with edges that are smooth apart from atomic scale steps. We find conductances 
quantized in integer multiples of 2e 2 jh and a plateau at ~ 0.5 x 2e 2 /h as in recent experiments 
[Tombros et at, Nature Physics 7, 697 (2011)]. We demonstrate that, surprisingly, conductances 
quantized in integer multiples of 2e 2 jh occur even for strongly non-adiabatic electron backscattering 
at the stepped edges that lowers the conductance by one or more conductance quanta below the 
adiabatic value. We also show that conductance plateaus near 0.5 x 2e 2 /h can occur as a result of 
electron backscattering at stepped edges even in the absence of electron-electron interactions. 

PACS numbers: 72.10.Fk,73.23.Ad,81.05.Uw 



Electrical conductances quantized in integer multiples 
of the fundamental quantum 2e 2 jh are the hallmark of 
ballistic quantum transport in nanostructures such as 
semiconductor quantum point contacts ji^ gold atomic 
wires r- and carbon nanotubes4 These quantized con- 
ductances are explained theoretically in terms of the 
Landauer theory of transport^ However, in the case 
of graphene nanostructures, quantum transport calcu- 
lations have shown the conductance quantization to be 
easily destroyed by disorder— ~— that is ubiquitous in 
these systems or by abrupt bends in the quantum wire 
geometry.— Accordingly, there have been only a few 



report: 



16-18 



of conductance quantization being observed 



experimentally in graphene nanostructures: Lin et al, le 
and Lian et alP demonstrated conductance quantiza- 
tion experimentally in graphene nanoribbons. However, 
the conductance steps that they observed were a few or- 
ders of magnitude smaller than the ballistic conductance 
quantum 2e 2 /h. This phenomeno n 16 i 17 has been ex- 
plained theoreticall y 1 ^ 20 as arising from strong electron 
backscattering at the edges of the electronic subbands of 
the ribbons due to the presence of random defects. More 
recently, Tombros et al^S- have reported the experimental 
observation of conductance quantization in integer multi- 
ples of 2e 2 /h, as well as a fractional conductance plateau 
at ~ 0.6 x 2e 2 /h, in a graphene nanoconstriction (GNC) 
at zero magnetic field. To minimize the effects of disor- 
der on transport in their device Tombros et alr^ stud- 
ied a short suspended GNC whose width was similar to 
its length and estimated to be ~ 300 nm. Their sam- 
ple was annealed by Joule heating which resulted in the 
constriction being formed with curved boundaries that 
were smooth on the mesoscopic length scale of ~ 100 nm. 
The atomic-scale structure of the boundaries was not de- 
termined experimentally, however, the curvature of the 
constriction's boundaries implies the presence of large 
numbers of atomic-scale steps (and possibly also other 
defects) along the boundaries. In this respect the GNC 
of Tombros et al^- differs from the well known semi- 
conductor quantum point contacts (SQPCs)^ where the 
transverse electron confinement is achieved electrostati- 



cally and thus the constriction boundaries are effectively 
smooth on the atomic scale as well as on the much larger 
(submicrometer) length scale of the overall dimensions of 
the constriction. In the limit of extremely slow spatial 
variation of the confining potential, electrons are adia- 
batically transmitted through the SQPC or adiabatically 
reflected. As was pointed out by Glazman et al^ such 
adiabatic transport results in quantized conductances; 
each electronic subband that is adiabatically transmitted 
through the narrowest part of the SQPC at the Fermi en- 
ergy contributes a quantum 2e 2 jh to the measured con- 
ductance. If the confining potential of the SQPC varies 
smoothly but not adiabatically, conductance quantiza- 
tion may still occur, each electronic subband transmitted 
through the narrowest part of the SQPC at the Fermi 
energy again contributing a quantum 2e 2 /h to the total 
conductance. 22 However, to date there have been no the- 
oretical studies of conductance quantization in constric- 
tions with boundaries exhibiting large-scale smoothness 
but atomic-scale steps, as in the GNC of Tombros et al*& 
For this reason a definitive understanding of the conduc- 
tance quantization in integer multiples of 2e 2 jh observed 
by Tombros et al^ has been lacking. Furthermore it has 
also been unclear whether the conductance plateau ob- 
served by Tombros et alJ£ at ~ 0.6 x 2e 2 /h was the result 
of electron-electron interactions (as is widely believed of 
the 0.7 x 2e 2 /h plateau in SQPCs^i) or whether it can 
be accounted for instead by strong single-electron scat- 
tering at steps in the constriction's boundaries. 

In this paper we report the results of quantum trans- 
port calculations that address these issues. We con- 
sider a non-interacting electron tight-binding model of 
graphene constrictions having similar dimensions to the 
GNC of Tombros et alJ£ and having boundaries that 
are smooth on the length scale of the constriction but 
with large numbers of steps on the atomic scale. We 
show that this model exhibits integer and fractional con- 
ductance plateaus similar to those that were observed 
experimentally. 18 Our results depend qualitatively on 
both the width of the constriction and its orientation. 
For the armchair orientation, the calculated integer quan- 
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FIG. 1: (color online) (a) Conductance vs. Fermi energy for 
GNCs with constriction widths W=150, 200, 250 nm. The 
constriction shape follows the cosine function and is smooth 
apart from steps on the atomic scale. Host and edge orienta- 
tion of the semi-infinite leads is armchair. Dotted lines show 
conductances of uniform armchair ribbons of the same widths 
as the narrowest parts of the GNCs whose conductances are 
plotted in the same colors. The black solid line shows the con- 
ductance of a 300 nm wide ribbon, whose width equals that of 
the leads W lead . (b) Conductance G vs. its energy derivative 
dG/dE. Dips in dG/dE indicate conductance plateaus. In 
the dashed rectangle dG/dE is shown for W — 200 nm only. 
Temperature T — 0. t — 2.7 eV. The subband spacing is an 
order of magnitude larger than ksT even at 4.2K as in Reflla. 



where e, is the on-site energy and tij =t = 2.7 eV is the 
matrix element between nearest-neighbor atoms. This 
Hamiltonian is known to describe the tt band disper- 
sion of graphene well at low energies^ Spin and elec- 
tron interaction effects are outside of the scope of our 
study. The nanoconstriction and any random edge dis- 
order and bulk vacancies that are present are introduced 
by removing carbon atoms from a uniform ribbon and 
setting appropriate hopping elements tij to zero. It is 
assumed that atoms at the edges are always attached to 
two other carbon atoms and all dangling bonds are passi- 
vated by a neutral chemical ligand, such as hydrogen, so 
that the bonding between the carbon atoms at the edge 
and around vacancies is similar to that in bulk graphene. 
Random bulk and edge disorder (when present) are char- 
acterized by the probability of the carbon atoms being 
removed, p and p e , respectively, p is normalized rela- 
tive to the whole sample, while p e is defined relative to an 
edge only. The long-range potential due to charged im- 
purities is approximated by a Gaussian form^iii of range 
d: Ci = J2 r Vbexp(— \ri — rrj| /d 2 ), where both the am- 
plitude Vo and coordinate tq are generated randomly. 

In the linear response regime the conductance of the 
GNC is given by the Landauer formula 13 
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Tji is the transmission coefficient from subband i in the 
left lead to the subband j in the right lead, at the Fermi 
energy. Tji is calculated by the recursive Green's function 
method, see Ref. [26|for details. The average conductance 
(G) for samples with random disorder was calculated by 
averaging over an ensemble of samples with different real- 
izations of the disorder. For the results presented below, 



tized conductances of the constrictions have smaller in- 
teger values than those of uniform armchair graphene 
ribbons with the same width as the narrowest part of 
the constriction. This differs qualitatively from the well 
known behavior of SQPCs where the adiabatic and non- 
adiabatic quantized conductance values are equal to those 
of a uniform quantum wire whose width equals that of the 
narrowest part of the constrictionj 21 ' 22 We find plateaus 
with conductance values ~ 0.5 x 2e 2 /h as well as the in- 
teger plateaus. For the zigzag orientation the calculated 
integer quantized conductances of the constrictions are 
either the same as or lower than those of uniform rib- 
bons of the same width as the constriction. We also find 
integer and fractional quantized conductances for con- 
strictions whose narrowest parts are neither zigzag nor 
armchair. 

We describe GNCs by the standard tight-binding 
Hamiltonian on a honeycomb lattice, 
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n\ FIG. 2: (color online) The same as Fig. [T]but for the zigzag 
configuration of host and edges in the leads. 



averaging was carried out over ten realization for each 
disorder type. 

To investigate the transport properties of GNC's 
we chose geometries similar to those studied 
experimentally. — The shape of constriction was modeled 
by a cosine function so that its edges were smooth apart 
from atomic scale steps. The width of narrowest part 
of the GNC was varied in the range W = 150. ..250 
nm. The GNC was attached at its two ends to semiin- 
finite leads represented by ideal nanoribbons of width 
W lead =300 nm. This guarantees that for any W the 
leads supply more states for propagation than can pass 
through the narrowest part of the constriction. The 
region of the constriction itself in our tight-binding 
quantum transport calculations included up to ~ 1.5 
million carbon atoms. In our modeling of the effects of 
random disorder we assumed it to be present only in a 
finite region of width 300 nm and length L = 300 nm; 
the semi-infinite leads were free from disorder. 

The calculated conductances of GNCs with different 
constriction widths W are shown in Figure QJa) f° r the 
armchair orientation of the graphene host and edges of 
the ideal leads. Note, however, that the edge orienta- 
tion along most of the constriction itself is neither arm- 
chair nor zigzag; see the outset in Figure [T] The con- 
ductance shows faint quantization steps in integer mul- 
tiples of 2e 2 /h, similar to those observed experimentally 
by Tombros et alM- For better visualization we plot the 
energy derivative of the conductance dG/dE in Fig. [TJb) . 
Here a dip in dG/dE indicates a plateau in the conduc- 
tance. The prominent dips in dG/dE in Fig. HJb) clus- 
ter around conductance values that are integer multi- 
ples of 2e 2 //i, including both odd and even integer multi- 
ples. The conductance of the GNC decreases as constric- 
tion becomes narrower, a feature expected theoretically 
and observed^ in conventional semiconductor quantum 
point contacts: As the constriction width shrinks the 
number of propagating states for a given Fermi energy de- 
creases. Interestingly, although the conductance plateaus 
occur near integer multiples of 2e 2 /h, in each case the in- 
teger has a smaller value than that for the ideal infinite 
ribbon of uniform width whose width equals the width 
W of the narrowest part of the constriction, for the same 
electron Fermi energy, calculated with the same tight 
binding approach. This can be seen by comparing the 
conductances of the GNC's in Fig. HJa) with those of the 
corresponding uniform ideal ribbon o ' 28 that are shown 
as the dotted lines of the same color in Fig. QJa). We 
also found no correlation between the calculated GNC 
conductances and the semiconductor/metallic property 
of uniform ideal armchair ribbons. These findings show 
that the conductance quantization that we find for the 
armchair oriented host and leads is not due to adiabatic 
transmission of individual eigenmodes of the ideal leads 
through the constriction but that additional scattering 
along the constriction edges plays an important role. 

For W = 200 nm, we find an additional conductance 
step at ~ 0.5 x 2e 2 /h; see Fig. [TJb) and the inset of Fig. 
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FIG. 3: (color online) Conductance vs. Fermi energy for GNC 
of width W=200 nm. The red open squares show the con- 
ductance for a constriction whose top and bottom parts are 
shifted by 80 nm relative to each other. Black filled squares 
show for comparison the conductance for the corresponding 
symmetric constriction of width W = 200 nm, as infjja). 



IJa). This agrees with the experimental findings in Ref. 
la . Whether or not this feature is present in the results of 
our quantum transport calculations depends on the width 
W of the constriction; note that experiments have been 
reported for only a single sampled However, as can be 
seen in Fig. [T]this fractional plateau coexists with integer 
conductance plateaus at higher electron Fermi energies 
for the same GNC and occurs for both electron and hole 
transport (not shown), as in the experimental dataJ^ 

Our results for GNCs with the zigzag orientatio n 29 i 30 
of the graphene host and edges of the leads are shown in 
Fig. [5J Comparison of Fig. Q]and Fig. [5] reveals signifi- 
cant differences between quantized conductance plateaus 
in GNCs with the host and leads in the zigzag and arm- 
chair orientations: For the zigzag orientation the quan- 
tized conductance plateaus are more pronounced than 
for the armchair case. Also for the zigzag case the calcu- 
lated values of the quantized conductances of the GNCs 
in many (but not all) cases are close to the values of the 
quantized conductances of ideal uniform zigzag ribbons 
having the same width as the narrowest part of the GNC 
and the same electron Fermi energy. By contrast, as 
we have already mentioned, all of the integer quantized 
GNC conductances for the armchair case are smaller than 
those of the corresponding uniform ideal ribbons by inte- 
ger multiples of 2e 2 /h. Thus in many cases non-adiabatic 
electron backscattering is much weaker for GNCs in the 
zigzag orientation than for those in the armchair orien- 
tation. This difference may be attributed to the current 
densities being much lower near zigzag graphene edges 
than near armchair edges so that the conductances are 
less affected by edge imperfections for zigzag ribbons. 31 

The open red squares in Fig. [3] show the calculated 
conductance of an asymmetric GNC with the armchair 
orientation of the host and leads. As shown in the inset, 
the geometry in this case is similar to the W = 200nm 
armchair-oriented constriction in Fig. Q] except that the 
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FIG. 4: (color online) (a) Effect of different disorder types 
on the conductance of armchair-oriented GNCs. (b) Conduc- 
tance averaged over 10 realization of disorder. Constriction 
width IP=200 nm. Red line with filled circles is for edge dis- 
order with p e = 0.2. Blue line with filled squares is for bulk 
vacancy disorder with p b — 10~ 5 . Green line with rhombuses 
is for long ranged potentials due to charged impurities with 
effective parameters \V\ < 0.2£, p = 5 x 10 15 m" 2 , d = 10a. 
The black solid line shows the conductance for an ideal, uni- 
form ribbon 300 nm wide and is given as a reference. 



upper and lower regions where the carbon atoms have 
been removed are now offset from each other laterally by 
80 nm. Thus the edges of the narrowest part of the con- 
striction have neither the armchair nor the zigzag orien- 
tation. We find that the electron backscattering is some- 
what stronger (the conductance lower) in this case than 
for the symmetric W = 200nm armchair-oriented con- 
striction in Fig. [TJ the calculated conductance for the 
latter is replotted as the solid black squares in Fig. [3] 
for comparison. However the first few quantized conduc- 
tance plateaus (as well as the plateau at ~ 0.5 x 2e 2 /h) 
are still clearly visible for the asymmetric GNC. 

The effects of disorder of different types are shown in 
Fig. [4] As a test system we chose a GNC of width 
W = 200 nm having the armchair orientation. The effect 
of disorder on the conductance of the GNC is similar to 
that for graphene nanoribbons. 19 However, the conduc- 
tance quantization is strongly degraded for every disorder 



type including bulk vacancies. This may be attributed to 
the varying width of the GNC along the transport direc- 
tion that precludes the existence of well-defined subband 
edges for the whole structure. We find each type of disor- 
der to suppress the conductance and to result in universal 
conductance fluctuations ! 19 ' 32 

In conclusion, we have carried out million-atom elec- 
tronic quantum transport calculations for graphene 
nanoconstrictions with boundaries that are smooth ex- 
cept for steps on the atomic scale and have dimen- 
sions similar to those of the graphene nanoconstrictions 
that have been found to exhibit conductances quantized 
in integer multiples of 2e 2 /h in recent experiments. 18 
Our results demonstrate quantized conductances simi- 
lar to those observed experimentally 18 in a tight binding 
model with non-interacting electrons. We find conduc- 
tances quantized in integer multiples of 2e 2 /h to occur in 
graphene nanoconstrictions even in the presence of strong 
electron backscattering at the stepped constriction edges 
that depresses the quantized conductance values by one 
or more 2e 2 /h conductance quanta below the quantized 
conductance values of uniform graphene ribbons with the 
same width and electron Fermi energy as those of the nar- 
rowest part of the constriction. This integer conductance 
quantization in the presence of such strong backscattering 
has no known analog in either adiabatic or non-adiabatic 
semiconductor quantum point contacts. It may explain 
why, based on their transport measurements, Tombros 
et al. 18 estimated the width of their GNC to be smaller 
(200-275 vs. 300 nm) at zero magnetic field than at 
higher magnetic fields where electron backscattering at 
the edges of the constriction is reduced. 3 — We also find 
that conductance plateaus at ~ 0.5 x 2e 2 /h need not be 
the result of electron-electron interactions in these sys- 
tems but can result instead from non-adiabatic backscat- 
tering of electrons at atomically stepped constriction 
boundaries. However, the plateau observed experimen- 
tally at ~ 0.6 x 2e 2 /h by Tombros et al~& resembles the 
plateau that is seen at ~ 0.7 x 2e 2 /h in SQPCs and is 
attributed to electron-electron interaction s 23 ' 24 in part 
because in SQPCs the potentials are smooth and there 
is no analog of the atomic steps present at the edges of 
GNCs. Therefore further experimental studies are re- 
quired to clarify whether electron-electron interactions 
or boundary scattering are primarily responsible for the 
fractional plateau observed by Tombros et alJ£ in the 
GNC. Our results (see Fig. [1]) suggest that systematic ex- 
perimental studies of GNCs having differing widths may 
answer this question. Our quantum transport calcula- 
tions also show random defects to strongly degrade the 
conductance quantization in graphene nanoconstrictions. 

This work was supported by NSERC, CIFAR, Com- 
pute Canada and WestGrid. 
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